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Abstract. We present a general method for the recon- 
struction of the one-point Probability Distribution Func- 
tion of the local aperture mass in weak lensing maps. Ex- 
act results, that neglect the lens-lens coupling and depar- 
ture form the Born approximation, are derived for both 
the quasilincar regime at leading order and the strongly 
nonlinear regime assuming the tree hierarchical model is 
valid. We describe in details the projection effects on the 
properties of the PDF and the associated generating func- 
tions. In particular, we show how the generic features 
which are common to both the quasilinear and nonlin- 
ear regimes lead to two exponential tails for P(M ap ). We 
briefly investigate the dependence of the PDF with cos- 
mology and with the shape of the angular filter. Our pre- 
dictions are seen to agree reasonably well with the results 
of numerical simulations and should be able to serve as 
foundations for alternative methods to measure the cos- 
mological parameters that take advantage of the full shape 
of the PDF. 
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1. Introduction 

Recent reports of cosmic shear detection (van Waerbeke 
et al. 2000a, Bacon, Refregier & Ellis 2000, Wittman et 
al. 2000, Kaiser, Wilson & Luppino 2000) have underlined 
the interest that such observations can have for exploring 
the large-scale structures of the Universe. Previous pa- 
pers have stressed that not only it could be possible to 
measure the projected power spectrum (Blandford et al. 
1991, Miralda-Escude 1991, Kaiser 1992) of the matter 
field, but also that non-linear effects could be significant 
and betray the value of the density parameter of the Uni- 
verse. More specifically, Bernardeau, van Waerbeke and 
Mellier (1997) have shown that the skewness, third order 
moment of the local convergence field, when properly ex- 
pressed in terms of the second moment can be a probe of 
the density parameter independently of the amplitude of 
the density fluctuations. This result can be extended to 



higher order moments, to the nonlinear regime (Jain & 
Seljak 1997, Hui 1999, Munshi & Coles 2000, Munshi & 
Jain 1999b) and the whole shape of the one-point PDF 
(Valageas 2000a,b; Munshi & Jain 1999a,b). 

In case of weak lensing surveys it appears however that 
it is more convenient to consider the so called aperture 
mass statistics that corresponds to filtered convergence 
fields with a compensated filter, that is with a filter of 
zero spatial average (Kaiser et al. 1994, Schneider 1996). 
Indeed, it is possible to relate the local aperture mass to 
the observed shear field only, whereas, in contrast, conver- 
gence maps require the resolution of a non-local inversion 
problem and are only obtained to a mass sheet degener- 
acy. The aperture mass statistics have proved valuable in 
particular for cosmic variance related issues (Schneider et 
al. 1998). Thus, in this article we present a method to 
compute the one-point PDF of the aperture mass, both 
for the quasilinear and strongly non-linear regimes. In the 
case of the quasilinear regime we can use rigorous pertur- 
bative methods while in the highly non-linear regime we 
have to use a specific hierarchical tree model (which has 
been seen to agree reasonably well with numerical simula- 
tions). Although the details of the calculations are specific 
to each case we point out the general pattern common to 
both regimes which is brought about by the projection 
effects. In particular, our methods are quite general and 
actually apply to any filters, though we are restricted to 
axisymmetric filters for the quasi-linear regime. Our re- 
sults for the non-linear regime, where there is not such a 
restriction, can also be extended to multivariate statistics 
(p-point PDFs). 

In Sect. 2 we recall the definitions of the local conver- 
gence and aperture mass and how they are related to the 
cosmic 3D density fluctuations. In particular we present 
the shape of the compensated filters that we use for the 
explicit computations we present in the following. In Sect. 
3 we describe the relationship between the PDF and the 
cumulant generating function of the 3D density field and 
we show how this extends to the projected density. The 
details of the calculations are presented in Sect. 4, for the 
quasilinear theory, and in Sect. 5 for the nonlinear theory. 
Numerical results are presented in Sect. 6. 
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2. The convergence and aperture mass fields 

In weak lensing observations, background galaxy defor- 
mations can be used to reconstruct the local gravitational 
convergence field. We recall here how the local conver- 
gence is related to the line-of-sight cosmic density fluctu- 
ations. As a photon travels from a distant source towards 
the observer its trajectory is perturbed by density fluctu- 
ations close to the line-of-sight. This leads to an appar- 
ent displacement of the source and to a distortion of the 
image. In particular, the convergence k magnifies (or de- 
magnifies) the source as the cross section of the beam is 
decreased (or increased). One can show (Kaiser 1998) that 
the convergence along a given line-of-sight is, 



0.6 



<m w(n,n s ) 5(n) 



(i) 



when lens-lens couplings and departure from the Born ap- 
proximation are neglected (e.g., Bernardeau et al. 1997). 
This equation states that the local convergence is obtained 
by an integral over the line-of-sight of the local density 
contrast. The integration variable is the radial distance, 
1Z, (and 1Z S corresponds to the distance of the source) 
such that 



dTZ = 



cdz/ Hq 



v^ A + (i - a„ - n A )(i + z) 2 + On(l + z? 

while the angular distance T> is defined by, 



V(z) 



c/H 



sinh ( \J\ — Q m — f^A 



HnTZ 



V 7 ! - — 

Then, the weight w(TZ,TZ s ) used in (|l|) is given by: 
30 m Hi V(K)V(K S - U) 



uj(K,K s ) = ^^ 



(1+*) 



(2) 



(3) 



(4) 



where z corresponds to the radial distance 1Z. Thus the 
convergence k can be expressed in a very simple fashion as 
a function of the density field. We can note from (jl|) that 
there is a minimum value K m i n (z s ) for the convergence 
of a source located at redshift z s , which corresponds to 
an "empty" beam between the source and the observer 
(5 = — 1 everywhere along the line-of-sight): 

f TZs 

K min = - dnw{n,n s ) (5) 

Jo 

In practice, rather than the convergence k it can be more 
convenient (Schneider 1996) to consider the aperture mass 
M ap . It corresponds to a geometrical average of the local 
convergence with a window of vanishing average, 



M„. 



d 2 0' U(&') k{0' - ■&) 



(6) 



where is the local convergence at the angular position 
# and the window function U is such that 



d 2 dU{$) = 0. 



(7) 
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Fig. 1. Shape of the filter functions we use in Fourier 
space. The solid line corresponds to the shape proposed by 
Schneider, Eq. (p^), multiplied by 1.459 and the dashed 
line corresponds to the compensated filter we introduce in 
this paper, Eq. ([IT 



In this case, M ap has the interesting property that it can 
be expressed as a function of the tangential component "ft 
of the shear (Kaiser et al. 1994; Schneider 1996) so that 
it is not in principle necessary to build local shear maps 
to get local aperture mass maps. More precisely we can 
write, 



M ap (tf) = / dVQO?')7t(0-0') 



with 



V 2 



di?' d'U {■&'). 



(8) 



(9) 



Nonetheless considering such a class of filters is interest- 
ing because convergence maps are always reconstructed to 
a mass sheet degeneracy only. Therefore, to some extent, 
any statistical quantities that can be measured in con- 
vergence mass maps correspond to smoothed quantities 
with compensated filters. In the following, for convenience 
rather than due to intrinsic limitation of the method, we 
consider filters that are defined on a compact support. 

2.1. Choice of filter 

It is convenient to write the filter function in terms of 
reduce variable, d/0 where 6 is the filter scale, 

li(0/0) 



(10) 



so that the evolution with 9 of the properties of M ap only 
depends on the behavior of the density field seen on dif- 
ferent scales (while the shape and the normalization of 
the angular filter Z7(i?) remains constant). In the follow- 
ing we shall use two different filters, which satisfy (|l0|). 
One, which we note us, has been explicitly proposed by 
Schneider (1996), 



us(x) = - (1 - x 2 ) 

TT 



for x < 1 



(11) 
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and ug(x) = otherwise. The Fourier form of this filter 
corresponds to, 



Ws(l) = 27r / Ax x us{x) Jq(Ix) 



24 J 4 (0 
I 2 ' 



The other, which we note ubv, corresponds to a simpler 
compensated filter that can be built from two concentric 
discs. It is built from the difference of the average conver- 
gence in a disc 1/2 and the average convergence in a disc 
unity: 



u BV (x) 



AH{2x) H{x) 



(13) 



where H (x) is the characteristic function of a disc unity. In 
Fourier space it is simply related to the Fourier transform 
of a disc of radius unity, 



Wi(l) 



2Ji(0 



and reads, 



W BV (l) = W 1 (l/2)-W 1 (l) 



4Ji(Z/2) 2Ji(Z) 



I 



I 



(14) 



(15) 



The ratio of the 2 radii has been chosen so that the 2 
filters are close enough, as shown in Fig. 1. Note that 
the normalizations are somewhat arbitrarily. In the plot 
they have been chosen to give the same amplitude for the 
aperture mass fluctuations in case of a power law spectrum 
with index n = —1.5. This is obtained by multiplying the 
expression @) by 1.459. Fig. 1 shows actually that the 
two filters are very close to each other. In particular they 
have their maximum at the same k scale, and except for 
the large k oscillations they exhibit a similar behavior. In 
the following we will take the freedom to use either one or 
the other for convenience. 



3. One-point PDF construction 

In this section we succinctly review the theory for the 
construction of the one-point PDF statistical quantities. 
In particular we recall the mathematical relationship be- 
tween the one-point PDF and the moment and cumulant 
generating functions. 

3.1. General formalism 

In general one can define xW as the generating function 
of the cumulants of a given local random variable 5, 



x (A) = X> P >c 



P =i 



(16) 



It is given through a Laplace transform of the one-point 
PDF of the local density contrast 5, 



e x(A) =y dSe xs P(S). 



(17) 



For hierarchical models, that is when (5 P ) C ~ (5 2 ) 
convenient to define tp(y) as 



it is 



(12) <p{v) 



E 



(S p )c (-y) p 



(<5 2 )p~ 1 pi 



-(s 2 ) x (-y/(S 2 ))- (is) 



Then the one-point PDF of 5 is then given by the inverse 
Laplace transform (see Balian & Schaeffer 1989), 



P(S)dS = d<5 / exp 
J 2ttig z 



v(y) , y$ 



(19) 



of the moment generating function. Here a = (S 2 ) 1 ^ 2 is 
the r.m.s. density fluctuation. 

3.2. The projection effects 

The relation ([!]) states that the local convergence can 
be viewed as the superposition of independent layers of 
cosmic matter field. The direct calculation of the one- 
point PDF of such a sum would involve an infinite num- 
ber of convolution products which makes it intractable. 
It is more convenient to consider the cumulant generat- 
ing functions which simply add when different layers are 
superposed (because Laplace transforms change convolu- 
tions into ordinary products). 

The projections effects for statistical properties of the 
local convergence have already been considered in previous 
papers. It has been shown in particular how the moments 
of the projected density can be related to the ones of the 
3D field in both hierarchical models corresponding to the 
non-linear regime (Toth, Hollosi and Szalay 1989) and in 
the quasilinear regime (Bernardeau 1995). 

More recently it has been shown (Valageas 2000a, b; 
Munshi & Jain 1999a) that these results could be extended 
to the full PDF of the projected density. In particular the 
cumulant generating function of the projected density can 
be obtained by a simple line-of-sight average of the 3D 
cumulant generating function. 

It is convenient to define the normalized projected den- 
sity contrast <5 pro j. by: 



Jproj. 



dKF(K)S(ll), 



(20) 



where F(1Z) is the selection function for the projection 
effects as a function of the radial distance 1Z: 



F(K) 



w(K,K s ) 



When filtering effects are included we have, 
<W,0= / d 2 d'w e ($') <W(0'-0) 



(21) 



(22) 



where wg is a given window function at scale 6. A par- 
ticular case is provided by the normalized aperture mass 



apt 



M ap = 



M, 



ap 



= / &&u{#)& vtol X#'-'d). 



(23) 
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The cumulants of the projected density can be related 
to those of the 3D density fields. Formally they correspond 
to the ones of the field when it is filtered by a conical shape 
window. Thus, from (BOI) we obtain: 



(Voj.(tfi)---<W(tffO>c = 



1Z S P 



Hdfti F{Hi) 



(24) 



The computation of such quantities can be made in the 
small angle approximation. Such approximation is valid 
when the transverse distances T>\&i — 1?| are much smaller 
than the radial distances TZ. In this case the integral 
( p4j ) is dominated by configurations where TZi — TZj ~ 
T>i\di — &j\ ~ X5j | i?i — It permits to make the change 
of variables TZi —> fi with TZi = TZ\ + r{D\. Then, since 
the correlation length (beyond which the many-body cor- 
relation functions are negligible) is much smaller than the 
Hubble scale c/H{z) (where H(z) is the Hubble constant 
at redshift z) the integral over r, converges over a small 
distance of the order of \di — and the expression ( p4~| ) 
can be simplified in, 

(*proj.(*i)---W(*p)>o = / ^rd^i-Ff 

Jo 

/oo P 
] { dr, ; (*(Wi,Z>i*i)...<y(^i+riX>i,2?i* p )) .(25) 
-°°i=2 

Taking filtering effects into account leads to, 

J proj.,9/c ~~ 



X 



l — l 

/CO f 
Y[d n (SCRuVitit) . ..S{U X + riZ>i,2?i* p )) c .(26) 
-°° i=2 



Thus the projection effects reduce to 



pro j • 



dTZF p (TZ) {5 V T 



where S^g cyl is the filtered 3D density with a cylindrical 
filter of transverse size T>8 and depth L (which goes to 
infinity in (|25|)). 

In particular this result gives the expression of the vari- 
ance of the filtered projected density contrast, 



(Si 



dTZF 2 (TZ) (Slg !Cyh ) c L 



(28) 



This expression can be re-expressed in terms of the power 
spectrum P(k) (nonlinear power spectrum), defined in this 
paper with 



(5( Xl )5(x 2 )) = 



d 3 k 

(2^)3 



P(k)e 



ik(xi-x 2 ) 



Then 



(<^I>e,cyl.)c 



dfc|| 
"2V 



d 2 k ± , N 



(2k) 



2 sin(fc|| L) 
fen L 



1 2 



W 2 {V9k 1 _) 



(30) 



where is the component of k orthogonal to the radial 
direction and k\\ is the component along the line-of-sight. 
In the previous integral, fc|| ~ 1/L and k± ~ 1/(D0), so 
that when L is large k\\ is negligible compared to k± which 
leads to, 



(6 2 



T>8, cyl. 



d 2 kj 



P(k±)W 2 (D6k ± ) 



(31) 



where W is the Fourier shape of the 2D window function. 
This relation holds for the filtered projected density con- 



trast as well as the aperture mass, for which W in (31) is 
to be replaced by Ws or Wbv- 

The formal expression for the higher order moments 
can be simplified by taking advantage of the so-called scal- 
ing laws for the correlation functions. It is in particular 
natural to assume that, 



(S p )c oc (8 2 ) 



2\P-1 



(32) 



with a coefficient of proportionality, S p , that depends on 
both the power spectrum and filter shapes, but not on the 
power spectrum normalization. For power law spectrum it 
implies in particular that these coefficients do not depend 
on the filtering scale. In the coming sections we present in 
more details the origin of this scaling relation. It allows to 
define 



<p(y) 



-{S 2 ) X (-y/(5 2 }). 



(33) 



The equation ( p7| ) then relates the cumulant generating 
function ip(y) for the projected density to the one corre- 
sponding to cylindrical filtering effects, 



(27) <p pX oj.(y) 



dTZ 

MR) 



<Pcyi.\yF(K)MK)] 



with 

(S 2 

MTl) = V, 2 



P9,cyl. 



L 



(34) 



(35) 



(Voj.,.) 

which can be rewritten in terms of the matter fluctuation 
power spectrum, 

/ d 2 kP(k,z) W 2 (kVt 



JdTZ'F 2 (TZ') J d 2 kP(k,z')W 2 (kV ( 



(36) 



In this expression we have explicitly written the redshift 
dependence of the power spectrum. In case of a power law 
spectrum, 

k\" 



(29) P(k,z)=P (z) 



k 



(37) 
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it takes a much simpler form given by, 

P (z)V- n - 2 



JdW F 2 {W)P (z')V'- 



(38) 



The result (|3J) is the cornerstone of the calculations we 
present. It allows to relate the cumulant generating func- 
tion of projected quantities to the ones computed in much 
simpler geometries. 

The difficulty then resides in the computation of tp cy \, 
in the regimes we are interested in. Two limit cases are 
actually accessible to exact calculations. First, one is the 
quasilinear regime where one can take advantage of the 
special properties of the perturbative expansion in La- 
grangian space to build the cumulant generating function 
in Eulerian space. These kinds of properties had been used 
previously for 3D or 2D top-hat filters only. We show here 



how this method can be extended to the filter (15). Sec- 



ond, in the strongly nonlinear regime one can also do ex- 
act numerical calculations when one is assuming the high 
order correlation functions to follow a tree model. The 
derivations corresponding to these regimes are presented 
in the next two sections. 

4. The quasi-linear regime 

The reason rigorous calculations can be carried on in this 
regime is that the cumulant generating function for a 
cylindrical shape is exactly the one corresponding to the 
2D dynamics (Bcrnardeau 1995). In other words, in this 
case 



quasilincar/ \ 



quasilincar/ \ 

<P2D (y) 



(39) 



at leading order. In this regime the problem then reduces 
to the computation of cumulant generating functions for 
compensated filters for the 2D dynamics. The latter cal- 
culation is presented in the following paragraphs. This is 
a long and technical calculation that leads to the formulae 



The calculations we present follow what has been done 
for the top-hat window filters, with the complication in- 
troduced by the use of two such filters instead of one to 
build the compensated filter. 

4--1- 2D statistics in Lagrangian space 

The generating function for the compensated filter ( |l3|) 
can be built from the generating function of the joint den- 
sity PDF for two concentric cells of different radius. This 
quantity will be obtained in Eulerian space from the one 
in Lagrangian space through a Lagrangian-Eulerian map- 
ping. More precisely we consider the joint PDF of two 
reduced volumes defined as the comoving volume V c {t) 
occupied by some matter expressed in units of the volume 
V it occupied initially, 



Vc(t) 

v • 



(40) 



The Eulerian overdensity of this matter region will then 
be given by the inverse of v, 



P = ^h- 



(41) 



The calculation will be made in two steps. First we present 
the derivation of the cumulant generating function in La- 
grangian space, then the mapping from Lagrangian to Eu- 
lerian space. 

In a Lagrangian description v corresponds to the Ja- 
cobian of the transform from the initial coordinates q in 
Lagrangian space to the ones in real space x, 



v = J(q) 



<3x 
<9q 



(42) 



The construction of the volume PDF is then based on 
the geometrical properties of the Jacobian perturbative 
expansion. Its expansion with respect to the initial density 
fluctuations (in the rest of this subsection we consider 2D 
dynamics) reads 



J(q) = l + J«(q) 



(43) 



Each term of this expansion can be written in terms of 
the initial Fourier modes of the linear density field, 



JW = 



d 2 ki d 2 k 



2tt 
exp[i(kx 



2^(*)*(ki). 
. +k p )q] J p (k l5 . 



.6(k p ) x 



■ j k p ) 



(44) 



where D + (t) describes the time dependence of the linear 
growing mode. The central issue is the way the geometri- 
cal kernel J p (ki, . . . , k p ) behaves when geometrical effects 
are taken into account. At leading order in perturbation 
theory (that is when only "tree order" terms are taken 
into account) that amounts to compute terms of the form 



J dai . . . da p J p (ki, . . . , k p ) W\ki 



(45) 



where cti is the angle of the i th wave vector. There exists a 
central property, valid for top-hat filters only, which states 
that (Bernardeau 1995), 

J d ai ...da p J p (ki,...,kp)Wi|k 1 + ... + kp| = (46) 

Wi (fci) . . . W% (k p ) J da!... da p J p (k% , . . . , kp) , 

where W\ is defined in (|l4|). This result extends the one 
obtained in Bernardeau (1994) for the 3D dynamics. 

At leading order in Perturbation Theory any cumu- 
lant of the form (vfv^jc involves only products of such 
quantities. The sort of commutation rule given in the pre- 
vious equation implies that these cumulants can be com- 
puted without explicitly taking into account the filtering 
effects. It means that any cumulant of the form {v\ u|) c 
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can be built with a tree shape construction with two differ- 
ent kinds of end points. Formally the generating function 
of such cumulants 



p.q 



p\ q\ 



reads, 



x(Ai,A 2 ) = AiO(ri) + A 2 0(t2) - 

ynCS(n)-^r 2 0(r 2 ) 
n = A iCnCj(n) + A 2 J 12 Cj(t2) 

r% = Aiei 2 CS(n) + A 2 ? 22 CS(r 2 ), 



(47) 

(48) 

(49) 
(50) 



where £y are the second moment of the linear density 
contrasts between two cells of fixed Lagrangian radii 0i 
and 9j, 



ijVn 3 j 



d 2 k 

(2tt)2 



P(fc) W(fc^) ^(fc^-)- (51) 



They are, for the Lagrangian variables Vi, fixed parameters 
that depend only on the power spectrum shape (and on 
the set of cell radii chosen at the beginning). The function 
£/ is the generating function of the angular averages of 



0( T ) = £ ip ~' 3o = 1 



p=0 



P 



3p 



(27T) 



K- J dai...dap J p (ki,...,kp 



(52) 
(53) 



Note that £(t) — 1 = 1/Qj(t) — 1 describes the density 
contrast of a spherical density fluctuation of linear over- 
density r for the 2D dynamics. The exact form of £j is 
therefore known for any cosmological model. 

The relation ( f48| ) can actually be generalized to an 
arbitrary number of cells in a straightforward way, 



n 



(54) 
(55) 



The latter relation can be rewritten in an equivalent way 

as 



Cy' r j — Aj Cj( T i) 



(56) 



where Cy is the inverse matrix to The generating 
function % can then be written, 



n ^ n 



(57) 



It gives the generating function of the reduced volume 
generating function of an arbitrary number of concentric 
cells. Note that however the known geometrical properties 
of the Lagrangian expansion terms do not allow to extend 
these results to non-concentric cells. Thus, our method ac- 
tually apply to any filter which is axisymmetric (in general 
one would need an infinite number of cells but in practice 
numerical discretization always leads to a finite number 
of concentric shells). It is worth noting that the relation 
(H) gives, 



3 



3 



which in turns leads to, 
dx({\}) 



(58) 



(59) 



4-2. Saddle point approximation and leading order 
cumulant generating function 

In this subsection we explicit the formal relationship be- 
tween the generating function computed at leading order 
and the shape of the multivalued density probability dis- 
tribution function. 

The joint PDF is formally given by (this is an extension 
ofEq.[l|) 



P(vi, 



. dAi dA n 
' 2n\ 2-7T1 



(60) 



x exp 



x({Ai})-$> 



In case of a small variance, the expression of the joint den- 
sity is obtained through the saddle point approximation. 
The saddle point conditions read, 



(61) 



which gives implicitly the values of at the saddle point 
position in terms of Vi. Taking advantage of the property 
^9h, one gets, 



(62) 



It implies that with the saddle point position the expres- 
sion of the joint PDF is (not taking into account prefac- 
tors), 



fW.(«i, 



exp 



ij 



(63) 



with the mapping (62). This is exactly what one would 



expect for Gaussian initial conditions, Tj being the linearly 
extrapolated density contrasts at the chosen scales. 
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4-3. Lagrangian-Eulerian mapping 

To relate Lagrangian and Eulerian space, one uses the 
same trick as in Bernardeau (1994), that is, 



Pi 



Vi > 



1 1 
, ...,v n > — 

P01 POn , 

-Peu1.G°1 < P01, • • • ,Pn < POn)- 



(64) 



The leading order cumulant generating function can then 
be obtained by an identification of the exponential term 
when one uses the saddle point approximation. The vari- 
ables are however now changed in pi which are related to 
T; with 



C(n) = Pi , C(n) = l/Ofa). 



(65) 



Moreover the variables pi enter also the expression of the 
cell correlation coefficients £y since they are in Eq. ( |5l| ) 
computed for a fixed mass scale and not for a fixed Eule- 
rian space radius. As a result, the coefficient Cy expressed 
in terms of £,„ should be understood as function of the 
variable pi through 



1/2, 



7 ( V 2 
Z{Pi 



V 2 a 



(66) 
(67) 



where 9i are all kept fixed. 

The cumulant generating function in Eulerian space is 
obtained also with a saddle point approximation in the 
computation of 

exp [xEui({Ai})] = J dpi . . .dp„ ffeuiGoi, ■ ■ ■ ,pn) 



x exp ^2 A»Pi 

which leads to, 

n ^ n 

XEui({A t }) = X iPi ~ 2 Gi i n T 3 
with the stationary conditions, 



2d Pi ^ °« 



A; 



(68) 



(69) 



(70) 



where the partial derivatives should then be understood 
for fixed radius 6i and Xi . The relation ( |69| ) , together with 
the conditions (f70| ) gives the formal expression of the cu- 
mulant generating function in Eulerian space. 
The case we are interested in, 



M ap = pi - p 2 
if 

01 = 1/2, e 2 = i 



(71) 
(72) 




-0.05 



0.05 



0.15 



Fig. 2. Shape of the functions Ti(y) and Tz(y) for the two- 
cell compensated window function (y = —Act 2 ). The com- 
putations have been done for a power law spectrum for 
n = —1.5. The solution for r has been extended beyond 
the singularity to show explicitly that the singularities are 
due to double solution in y. 



corresponds to 2 cells, if M ap is built with the filter (|15[). 
Then the generating function for M ap is obtained with a 
peculiar choice for Aj, 



A — Ai — — A2. 

It is actually convenient to define, 

y = -\cr 2 

and 

<P2D(y) = VoylXv) = -X(A)<7 2 

where 1 is defined by 

a 2 = (M a 2 p >=Cii+?22-2Ci2- 



(73) 
(74) 
(75) 
(76) 
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With this choice of variable, tp(y) does not depend for- 
mally on the variance but only on y. To be more specific 
one has finally, 



fcylXv) = y({Tl)-y((T2)+J r (T 1 ,T 2 ) 



2[l-r(n,7a)] 



£ll( T l) &22O2) 



2r(ri,T 2 )riT2 

?il(n)?22(T2) 



(77) 



(78) 



where r = £12/ \/ Ch^22j £11 an d £22 are considered as 
function of n and T2 through the variables p\ and p 2 . the 
saddle point conditions then read, 



dT 



(79) 
(80) 



In this case the function can be numerically cal- 
culated. We have restricted our calculations to the case 
where the power spectrum follows a power law behavior 
with index n = —1.5. To do the numerical computations 
we also use a simplified expression for the 2D spherical 
collapse dynamics, 



C(t) = (1+t/«T k 
with 



vT3- 1 



1.30. 



(81) 



(82) 



The resulting function tp C y\. (y) is shown in Fig. g, together 
with the functions T\(y) and T%(y). 

4-4- Properties of the cumulant generating function 

These figures clearly show that the function tp(y) has two 
singularities on the real axis. This is to be compared to 
what is encountered for counts-in-cells statistics where 
only one singular point is expected. The numerical resolu- 
tions have been extended slightly beyond the singularities 
to show that they are due to the resolution of the implicit 
equations in r that have multiple solution in y. As a result 
the generic behavior near any of such singularity is 



r(y) ~ t.b-y.) 1 '*, 
^cyi.(y) ~ a s (y - y s ) 3/2 . 



(83) 
(84) 



This behavior directly induces exponential cutoffs for the 
shape of the density PDF (see Balian & Schaeffer 1989, 
Bernardeau & Schaeffer 1992). In case of compensated fil- 
ter, the fact that we obtain 2 singularities, induces two 
exponential cut-offs on both side of the PDF as it ap- 
pears clearly on the results presented in Sect. ^ (see also 
Valageas 2000c). 



U (t) 




Fig. 3. Geometrical representation of the equation flS6|). 
The function r(y) is given by the intersection of the line 
—r/y with the curve C,' eS (r). Singularities appear when 
both curves are tangential, as shown in the figure. 

In a phenomenological way the function <p cy \.{y) can 
be described by an effective vertex generating function 
Cofr.(r) so that, 

Vcyi.(y) = yCeflE.(-r) - ^yrCeff.(-r) (85) 

t = -vC S .(r). (86) 

Numerically the effective vertex generating function is well 
described by a fifth order polynomial, 



Ceff.(T) = -T- 



0.843179 t ■ 
6.3852 r 5 



5.7511 r 3 + 3.3669 r 4 



(87) 



which is regular, the expected singular behavior for ip(y) 
being induced by Eq. ( p6| ) (see Fig. ||). Note that (eff.( T ) 
can be viewed as a generating function of effective vertices. 
It implies that for instance the skewness of the 2D (or 
equivalently cylindrical) compensated filtered density is 



s: 



cyl. 



3 x 0.843179 « 2.52 



(88) 



for such power law spectrum shape. It is however impor- 
tant to have in mind that the shape of £ G ff . (r) as well as the 
skewness depend on the window function normalization. 
The calculations have been given here for the ubv filter 
defined in Eq. @. For the u s filter, Eq. @, 5g yL for in- 
stance would have been about 1.459 times larger, because 
of the normalization discrepancy between the two filters. 

The skewness of M ap is then related to this one through 
a simple projection factor, 

, cyl JdHF 3 [P (*)P-("+ 2 >] 2 



5'| >roj 



5': 



[J dKF 2 P Q (z)V-( n +V] 

This latter relation is actually valid in both the quasilinear 
and the nonlinear regime. 

5. The nonlinear regime 

In the nonlinear regime, there exist no derivations from 
first principles of the behavior of the high order correla- 
tion functions of the cosmic density field. However, hints 
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of its behavior can be found. The stable clustering ansatz 
gives indication on the expected amplitude of the high or- 
der correlation functions and how they scale with the two- 
point one. The hierarchical tree model, and more specif- 
ically the minimal tree model, allows to build a coherent 
set of high-order correlation functions. 

5.1. The stable- clustering ansatz 

The stable clustering ansatz (Peebles 1980) simply states 
that the high order correlation functions should compen- 
sate, in virialized objects, the expansion of the Universe. 
It gives not only the growth factor of the two-point corre- 
lation, but also a scaling relation between the high-order 
correlation functions. 

Expressed in terms of the coefficient S p -hence of the 
generating function tp(y) defined in (p"8|)- it means that 
they are independent of time and scale. As a consequence, 
the knowledge of the evolution of the power-spectrum 
P(k), or of the two-point correlation function £, is suffi- 
cient to obtain the full PDF of the local density contrast. 
This property has been checked in numerical simulations 
by several authors (Valageas et al. 2000; Bouchet et al. 
1991; Colombi et al. 1997; Munshi et al. 1999). In partic- 
ular, the statistics of the counts-in-cells measured in nu- 
merical simulations provide an estimate of the generating 
function tp(y) for 3D top-hat filters. 

More precisely, in the highly non-linear regime one con- 
siders the variable x defined by: 

l + S R 



(90) 



Then, using (Jl^), for sufficiently "large" density contrasts 
the PDF P(8r) can be written as (Balian & Schaeffer 



(91) 



P(Sr) = |j K*) 
when 

e»i, (i + ^)»e 



-u/(l—u>) 



(92) 



where the scaling function h(x) is the inverse Laplace 
transform of (p(y): 



/-Hog 7 



(93) 



In (^2|) the exponent u> comes from the behavior of (f(y) at 
large y. Indeed, from very general considerations (Balian 
& Schaeffer 1989) one expects the function ip(y) defined 
in ( |l8| ) to behave for 3D top-hat filtering as a power-law 
for large y: 

y +oo : tp(y) ~ a y x ~ u with < w < 1 , a>0 (94) 
and to display a singularity at a small negative value of y, 
!/->!/+ : tp(y) = -a. T(lu s ) (y - y s )~ u * (95) 



where we neglected less singular terms (note that this be- 
havior has indeed been observed in the quasilinear regime, 
Bernardeau 1992). Taking advantage of these assump- 
tions, one obtains (Balian & Schaeffer 1989), 



, / \ °(1 — W ) ,,9 

x < 1 : h{x) y ' x u ~ 2 

1 (u) 

x 1 : hix) ~ a s x^ B * e 



(96) 
(97) 



with x s — l/\y s \ - Hence, using ( |9l| ) we see that the density 
probability distribution P{6r) shows a power-law behav- 
ior from (1 + Sr) ~ £ ^ ' up to (1 + 5r) ~ x s £ with 
an exponential cutoff above x s ^. It implies in particular 
that the function h(x) measured in numerical simulations 
can give rise to constraints on the cumulant generating 
function from the inverse relation, 



1 



cy ) h(x) dx 



(98) 



Note that h{x) depends on the power-spectrum and, in 
the absence of a reliable theory for describing the nonlin- 
ear regime, it has to be obtained from numerical simula- 
tions. This is the case in particular for the evolution of 
the two-point correlation function, or equivalently of the 
power-spectrum. To this order we use the analytic for- 
mulae obtained by Peacock & Dodds (1996) from fits to 
N-body simulations. 

Note that the relation ( |9^ ) holds independently of the 
stable-clustering ansatz. However, if the latter is not re- 
alized the generating function (p(y) depends on time (and 
scale). Then, most of the results we obtain in the next 
sections still hold but one needs to take into account the 
evolution with redshift of tp(y). That would be necessary 
in particular if one wants to describe the transition from 
the quasilinear regime to the strongly nonlinear regime. In 
the following we assume that the stable-clustering ansatz 
is valid, so that (p(y) is time- independent. As mentioned 
above this is consistent with the results of numerical sim- 
ulations. 

5.2. Minimal tree-model 

If one is interested in the statistics of the top-hat filtered 
convergence, it is reasonable to assume that (Valageas 
2000a,b), 



^cyl.(2/) ~ <P3D(y)- 



(99) 



In case of the aperture mass statistics however the filtering 
scheme is too intricate (with both positive and negative 
weights) to make such an assumption, and in particular 
the resulting values for S p depend crucially on the geo- 
metrical dependences of the p-point correlation functions. 
We are thus forced to adopt a specific model for the cor- 
relation functions, and the one we adopt is obviously con- 
sistent with the stable-clustering ansatz. 
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A popular model for the p— point correlation functions 
in the non-linear regime is to consider a "tree-model" 
(Schaeffer 1984, Groth & Peebles 1977) where £ p is ex- 
pressed in terms of products of £2 as: 



^(r 1 ,...,r p )=^Q(«)^n^(r l ,r J ) 



(100) 



(a) 



p-1 



where (a) is a particular tree-topology connecting the p 
points without making any loop, Qp is a parameter asso- 
ciated with the order of the correlations and the topology 
involved, t a is a particular labeling of the topology (a) 
and the product is made over the (p — 1) links between 
the p points with two-body correlation functions. A pecu- 
liar case of the models described by ( |100| ) is the "minimal 
tree-model" (Bernardeau & Schaeffer 1992, 1999, Munshi, 
Coles & Melott 1999) where the weights Qp ' are given 
by: 



n - 

vertices of (a) 



(101) 



where v q is a constant weight associated to a vertex of the 
tree topology with q outgoing lines. Then, one can derive 
the generating function cp(y), defined in (|Ts|), or the co- 



efficients S p , from the parameters v p introduced in (101) 
which completely specify the behavior of the p— point cor- 
relation functions. 

In this case the cumulant generating function is given 
for 3D filtering by, 



X(A) = 



P =i 



d 3 ri ..d 3 r p w{r 1 )..w(r p ) f p (n, ., r p )(102) 



where w corresponds to the filter choice. In the case of 
the minimal tree-model, where the p— point correlation 



functions are defined by the coefficients v q from (|100|) and 



(101), it is possible to obtain a simple implicit expression 
for the function x(A) (see Bernardeau & Schaeffer 1992; 
Jannink & des Cloiseaux 1987): 



X(A) 



A 



J d 3 r w;(r) 



C[r(r)] 



r(r)C'[r(r)] 



(103) 



r(r) = \J d 3 r' w(r') £ 2 (r,r') C'[r(r')] (104) 
where the function £(t) is defined as the generating func- 



tion for the coefficient v„ 



C(r) 



E 

p=0 



(-iy 



pi 



V p T P with Vq = V\ = 1. 



(105) 



The function x(A) obviously depends on the choice of fil- 
ter through the function w. For a top-hat filter, it would 
simply be a characteristic function normalized in such a 
way that 



d 3 r w(r) = 1. 



(106) 



A simple "mean field" approximation which provides 
very good results in case of top-hat filter (Bernardeau & 
Schaeffer 1992) is to integrate r(r) over the volume V 
in the second line of the system (103) and then to ap- 
proximate r(r) by a constant r. This leads to the simple 
system: 



<P3£>(y) = y 



C(r) 



rC(r) 



r = -yC(r) 



(107) 
(108) 



Then, the singularity of tp(y), see (|95|), corresponds to the 
point where the |dy/dr| vanishes. Note that £(r) is regu- 
lar at this point and that the singularity is simply brought 
about by the form of the implicit system ( |l08| ) as observed 
in Bernardeau & Schaeffer (1992). Making the approxima- 
tion (107- 108 ) for both ^315 (y) and ip cy i.(y) leads to the 
approximation (^) which is thus natural for the minimal 
tree model. 

In the case of a compensated filter such a simple mean 
field approximation however cannot be done. It is in par- 
ticular due to the fact that the weights given to r then 
strongly depend on the radius distance. Before we go to 
this point we need first to take into account the projection 
effects. 

5.3. Projection effects for the minimal tree-model 

As noted in Toth et al. (1989) and analyzed in detail 
in Valageas (2000b), we know that the tree structure as- 
sumed for the 3D correlation functions is preserved (ex- 
cept for one final integration along the line-of-sight) for 
the projected density. Indeed, inserting (10C) in fl24| ) we 
obtain: 



(Wtfi)--Wtf P ))c = ]TQ( Q) 



E 

(a) ta 

p 



dlli F p 



/DO P 
H d ^ II>(Xa,X6) ( 109 ) 
-OO -_o „ 1 



p-i 

where we noted x a = (lZ a ,T>i-d a ). It can be noted that 
in the small angle approximation the weight applied to 
each diagram depends on their order p only and not on 
their geometrical decompositions. As a consequence the 
projected p— point correlation function can be written, 

(<Wj.(0i)---<W(0 P ))c = 



dlZ F p u p (0 1 ,...,0 p ;z), 



(110) 



where the two-dimensional p— point functions 
£jj,(i?i, 'dp; z) have the same tree-structure as the 
three-dimensional p— point correlation functions £ p , 

wp^i, . . . = Q P a) E II ^(tfa, z ) ( m ) 

(a) ta P-1 
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with: 



w 2 (tfi,tf 2 ;z) = J-^P(k,z) J o (*Z>|0i-0a|). (112) 

Here Jo is the Bessel function of order 0. For convenience 
we also note 072(2) the angular average of ct^C^i, #2; 2), 

57 2 («) = J d 2 'di& 2 'd2U{ti 1 )U{ti 2 )uj 2 {'d 1 ,ti 2 -z), (113) 

which, expressed in terms of the power spectrum gives, 




w 2 (z) 



d 2 k 

(2tt)2 



P(fc) W 2 (fcX>(9). 



(114) 



Thus, we see that the correlation functions of the pro- 
jected density <5 pr oj. itself do not show an exact tree- 
structure as the underlying 3D correlation functions £ p . 
Nevertheless, as seen in ( [L10[ ) they are given by one sim- 
ple integration along the line-of-sight of the 2D p— point 
functions u p which exhibit the same tree- structure as their 
3D counterparts £ p . This means that we can still use the 
techniques developed to deal with such tree-models. In 
particular, in the case of the minimal tree-model ( |101 ^ we 



0.4 0.6 
radius, d/6 

Fig. 4. Profile of the r function as a function of the an- 
gular radius. The two plots correspond to the two singu- 
larities, y s ~ for the solid line and y s + for the dashed line. 

As a consequence, in the nonlinear regime, the relation 
(p3) is to be used with, 



fcylXv) = y / d 2 tf U(#) 



C[r(0)] 



r(tf)C'[r(tf)] 



(117) 



will be able to take advantage of the resummation (102- 



104). 



r(*) = -„ / dV U{*) ^fy Z) " 



CW)] (118) 



Once again, it is interesting to note that for power- 
law spectra, P(k) oc k n , the angular and the redshift de- 
pendences of lo p can be factorized so that the correlation 
functions of the projected density S plo j. itself now exhibit 
a (new) tree-structure. Then, the 2-point function reads, 

(W(*i)W(0 2 )> = |*i-* 2 |- (n+a) 

roo j2i r p 

l n J o (0 / dK F 2 (K) p-("+ 2 ) (115) 

'0 



(2tt) 2 " " uy " J J k, 



where £(r) is the 3D vertex generating function. Note that 
this function depends on z through ol> 2 ($, 1?'; z) and u7 2 (;z). 
Note also that in y = the expansion of the generating 
function ip C y\.{y) is <p cy \.{y) = -J/ 2 / 2 + . . .. 

We have computed the resulting shape of the generat- 
ing function in such a model for various cases. For compar- 
ison with the previous quasi-linear case we assume here the 
power spectrum to follow a power law behavior with index 
n = —1.5. The vertex generating function is assumed to 
be given by 



while 



the 



high-order p— point functions C( T ) — (1 + T / K ) 



(119) 



(a) 



(<5proj.($l)--£pro.j.($j)))c follow the tree-structure (100) 
with the projected weights Q p ' 



'^ proj - ^ P \JdKF*P (z)T>-l<**)] 



Note that the relation Q p a ^ <-* Qp"p ro j depends on the 
slope n of the power-spectrum. On the other hand, if the 
initial tree-model for the 3D correlation functions is the 
minimal tree-model ( |l0l| ) we can see that the projected 
tree-structure ( |116| ) is not an exact minimal tree-model]^] 
which would be expressed in terms of a new generating 
function Cproj.^)- In other words cp pro j.(y) cannot be built 
from a tree structure whereas fcyi.(y) cai b an< l with the 
same vertex generating function £(r) defined in (105). 



p-i 



p-i 



(116) 



with hi ~ 0.5 (in the parameterization of £ we followed the 
traditional notation and used simple free parame- 

ter. It is not to be confused with the local convergence). 
In Fig. H we present typical profiles obtained for r(i?). We 
see that it is regular in §. In particular it does not exhibit 
discontinuities nor abnormal behavior near the singular 
values of y. We also found that the results we obtain are 
very robust regarding to the number of shells used to de- 
scribe the integral in with 2 cells only the description 
of tpcyi. is already very accurate. 



1 T his is due to the fact that the numerator in the r.h.s. of 
(11C) cannot be factorized in the form AB P_1 . A simple way 
to check that Qp Q p ro j cannot be written in terms of new pa- 
9,proj. as in (101) is to consider the "snake" topology 



The choice of the value of k in (119) relies a priori on 
numerical results. The EPT (Colombi et al. 1997) or more 
convincingly the HEPT (Scoccimarro & Frieman 1999) 
provide however a convenient frame which can be used 
to predict the value of k. This can be done for instance 
by identifying the predicted values for the skewness both 
from the form ( |119| ) and HEPT. Indeed in our model we 
have, 



rameters v. 



where Q p 



2 p-2 

v\ v 2 ■ 



S 3 (k) 



,(! + «) 



(120) 
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whereas, in HEPT, S3 is related to the initial power spec- 
trum index, 



->HEPT 



(n)=3 



4-2™ 
1 -(- 2"+ 1 ' 



which leads to, 



(l + 2 1+ ™ 



(121) 



(122) 



3 (1-2™) 

In the numerical applications presented in the following 
we will use this scheme. In particular, n = —1.5 leads to 
K = 0.88. On the other hand, at the angular scale 9 = 
4' which we consider below the local slope of the linear 
power-spectrum is n ~ —2.2 which leads to k ~ 0.6. 

The properties of ip C yh(y) we get in this regime are 
very similar to those obtained for the quasilinear regime. 
In particular we found that the function tp cy \. (y) exhibits 
2 singular points on the real axis. As for the quasilin- 
ear regime this behavior is due to the implicit equation 
in t and not to peculiar choice of the vertex generating 
function. In Table [l] we summarize the parameters that 
describe the singularities of <f C y\.{y) m different regimes. 
It appears, as expected, that the singularities are closer to 
the origin. This is to be expected since the nonlinearities 
contained in ip cy i, (y) are stronger in the nonlinear regime 
compared to the quasilinear regime. 

Similarly to the quasilinear regime it is also possible to 
define an effective vertex generating function from which 
Vcyi.(y) can be built and which reproduces its singular 
points, 



Ceff.(r) 



-t+ 1.4966 r 2 - 11.6982 r 3 
77.1899 r 5 . 



21.528 r 4 



(123) 



The resu lt given here has been obtained for k = 0.88 in 
Eq. (|ll|). 



6. Statistics of the aperture mass M. 



ap 



It now suffices to plug the numerical expressions we have 
obtained for ip cy \, (y) in rtjij) to get the shape of the M ap 
PDF. More precisely we have, 

dy 



P(M ap )dM a 



x exp 



I dAf ap 



2m<j 2 



(y) , 



ap 



(124) 



where a is the variance of the aperture mass. 

In (124) the integral over y has to be made in the com- 
plex plane. The integration path in the y plane is built in 
such a way that the argument of the exponential is always 
a real negative number thus avoiding oscillations (see Fig. 
|5|). Moreover, one must make sure that the integration 
path does not cross the branch cuts of tp cy \.(y). The sin- 
gularities of tfcylXv) induce non-analytic parts for ip plo j.(y) 
as well. They are located at positions, 



V > Vs 



and y < y%™' 



(125) 



0.02 



0.00 



-0.02 0.00 0.02 

real part, of y 



0.04 



Fig. 5. Integration paths in the y complex plane. The 
paths (dashed lines) are dynamically built so that the ar- 
gument of the exponential in Eq. (124) is always a negative 
real number. The thick half straight lines represent the lo- 
cations of the non-analytic parts of <^ pr oj.(j/) hi this plane. 
Of course, the paths must not cross these branch cuts. 



with 



j£± J =Wmax[F(7e)V(7e)] 



(126) 



where the maximum value of F(TZ)ip(lZ) is taken along the 
line-of-sight (and is indeed finite). As noticed in Valageas 
(2000a) the exponent ui s of the singularity of ip C yl.(y) 
(as defined in (9^)) leads to the exponent uo s — 1/2 for 
the projected generating function </5 pro j.(j/). However, in 
both the quasi-linear and highly non-linear regimes we 
have lo s = —3/2, see (84). In this shown in 

App. A, for y — > y^± 3 the singularity is of the form 

<^ P roj.(2/) ~ (y - yT± 3 ) 2 ln I y - yT± 3 \- 

The existence of these branch cuts is directly respon- 
sible for two exponential cut-offs in the shape of the PDF 



of M, 



ap ■ 



P(M ap ) - exp 



''min I -A^ap 



Vs± 



(127) 



It can be noted that the £1 dependence of K m i n will induce 
a strong dependence in the position of the exponen- 
tial cut-offs. The variation of re m i n with $1 is thus to be 
compared with the theoretical uncertainties on y s ±. 

6.1. The PDF shape 

In Figs. |6| we present the resulting shape of the one point 
PDF obtained for different cases. They have been obtained 
from the parameterization of tf cy \.{y) described in the pre- 
vious sections. In particular we assume a power law spec- 
trum with index n = —1.5. The variance adopted for the 
plots is a = 0.01. In these investigations we did not try to 
put a realistic source distribution, but we assume all the 
sources to be at redshift unity. However, all our results can 
be extended in a straightforward fashion for any redshift 
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Table 1. Values of S^ yL and singularity positions of <p C y\.{v) f° r Ubv filter. 

S^ y1 ' y a - y s + a s - a s + 

quasi-linear, n = -1.5 2.53 -0.092 0.159 1.20 0.74 

non-linear, k = 1.0, n = -1.5 4.19 -0.062 0.130 0.85 0.41 

non-linear, k = 0.88, n = -1.5 4.49 -0.057 0.121 0.83 0.40 

non-linear, k = 0.50, n = -1.5 6.32 -0.040 0.084 0.65 0.31 
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Fig. 6. Shape of the one-point PDF of the aperture mass 
from the quasi-linear regime model (top panel) and in the 
Nonlinear regime (bottom panel). The sources are at red- 
shift unity, the variance of M ap is 0.01, for the filter given 
by Eq. ITS the solid lines correspond to Einstein-de Sitter 
case, the dashed lines to flat universe with f2 m = 0.3. 

distribution of the sources (the latter is simply absorbed 
by a redefinition of the selection function F(1Z)). Note 
also that we present the PDF for the correctly normal- 
ized aperture mass M ap , that is without dividing the local 
convergence by K m i n - We can check in Fig. || that the tails 
of the PDF are stronger in the non-linear regime than in 
the quasi-linear regime independently of the variance a 
(which is the same in both plots). This is related to the 
smaller values of the singularities \y^+ 3 | in the non-linear 
case, as shown by the expression (127). Of course, this is 
due to the smaller value of the singularity \y s \ of the 3D 



Fig. 7. Same as previously with the nonlinear model but 
with cell radius ratio of 10 instead of 2. 

density field in the non-linear regime (in a similar fashion 
the coefficients S p are larger). 

In Fig. [j] we show that the positions of the cut-off de- 
pend crucially on the window shape. In particular it is 
clear that when the disc radius ratio is larger the PDF 
is more asymmetric and bears more resemblance with the 
k-PDF for a top-hat window function. Indeed, when the 
inner disk is much smaller than the outer radius the fluc- 
tuations of M ap = Ki — k,2 are dominated by those of 
the convergence K\ which corresponds to this small inner 
window while K2 which is governed by larger scales shows 
lower amplitude fluctuations. 

6.2. The f2 m dependence of the PDF 

It has been stressed in the literature (e.g., Bernardeau 
et al. 1997) that the non-Gaussian properties of the con- 
vergence maps are expected to exhibit a strong J7 m de- 
pendence. This is due in particular to the normalization 
factor. Such dependence is apparent in K m i n that depends 
crucially on the value of f2 m . This property naturally ex- 
tends to the shape of the one-point PDF. In particular it 
is important to have in mind that the f2 m dependence is 
negligible in ip cy \. (y) (for a fixed shape of the power spec- 
trum) . The Q m dependence is therefore entirely contained 
in the projection effect through the shape and amplitude 
of the efficiency function. 

In Fig. ^ we show how low values of r^ m amplify the 
non-Gaussian features contained in the PDF. Whether 
such a parameter can be constrained more efficiently with 



14 



F. Bernardeau & P. Valageas 



the PDF than with simply the local skewness is not yet 
clear. Such a study is however beyond the objective of this 
paper and is left for further works. 

6.3. Comparison with numerical simulations 

Finally, we compare our predictions for the PDF P(M ap ) 
of the aperture mass M ap with the results of N-body sim- 
ulations (Jain, Seljak & White 2000) in Fig. |. Note that 
for all these comparisons we exclusively use the us filter. 
We consider the cosmological models defined in Tab.||: a 
standard CDM (SCDM) and a rCDM scenario in a criti- 
cal density universe, a low-density open universe (OCDM) 
and a low-density flat universe with a non-zero cosmologi- 
cal constant (ACDM). Here T is the usual shape parameter 
of the power-spectrum. We use the fit given by Bardeen 
et al. (1986) for P(k). We only consider the weak lensing 
distortions which affect a source at redshift z s = 1, with 
angular window characteristic scale 9 = 4'. 

In the numerical calculations, we discretize the i n tegra l 
( |34| ) over redshift and we solve for the system ( 117 - |ll8| ) . 



That is we take into account the redshift dependence of the 
generating function (p cy \.(y). Moreover, we use the relation 
(122) to get the value of the parameter k, where for n we 
take the local slope of the linear-power spectrum at the 
wavenumber k = 2/ ' (V(z s /2)9). This corresponds to the 
Fourier modes which are probed by the filter of angular 
radius 9, see Fig. [p. For 9 — 4' we obtain n ~ —2.2 and 
k ~ 0.6. 

First, we can check in Fig. || that we recover the right 
trend for P(M ap ), with two asymmetric tails for large 
| M ap | . In particular, the exponential cutoff is stronger for 
negative values of the aperture mass than for positive val- 
ues. We can also note that the PDF is significantly differ- 
ent from a Gaussian as it shows a clear exponential cut- 
off, much smoother than the Gaussian falloff, especially 
for large positive M ap . On the other hand, we note that 
Reblinsky et al. (1999) obtained a good match to the tail 
of the PDF P(M ap ) using a description of the density 
field as a collection of virialized halos (Kruse & Schnei- 
der 1999). However, such a method is restricted to the far 
tail of the PDF (large positive M ap ) while our approach 
provides in principle a model for the full PDF P(M ap ). 
There seems to be a small discrepancy with the simula- 
tions for the rCDM scenario. It is not clear whether this 
is due to a limitation of HEPT or of our formulation. To 
clarify this problem one should test the statistics of the 3D 
density field and of the projected density in the same sim- 
ulation. However, this is beyond the scope of our paper. 
Nevertheless, the overall agreement appears to be quite 
reasonable. Note that the shape of the PDF is governed 
by only one parameter k, which is uniquely related to the 
local slope of the power-spectrum, independently of scale 
and of the cosmology. On the other hand, the inaccuracy 
of the numerical simulations in the tail of the PDF might 
be somewhat underestimated. 




Fig. 9. The relative difference [P(M ap ) - 

G(M ap )]/G(M ap ) of the aperture mass PDF P(M ap ) 
with respect to the Gaussian G{M ap ). The variance of 
the Gaussian is taken from the simulations. As in Fig. || 
we consider a source at redshift z s = 1 with the angular 
window 9 — 4'. We display our results for the rCDM 
(solid line) and OCDM (dashed line) cosmologies. The 
dotted lines correspond to the same cosmologies with 
K = 0.88 (see main text). The points show the results of 
N-body simulations from Jain, Seljak & White (2000). 

In order to see more clearly the difference of the aper- 
ture mass PDF P(M ap ) with respect to the Gaussian we 
display in Fig. || the relative difference P(M ap )/G(M ap ) — 
1. Here G(M ap ) is the Gaussian with the same variance as 
the numerical simulations. We can check that we recover 
a reasonable agreement with the numerical results, since 
Fig. [)] is directly related to Fig. ||. To get an estimate of 
the sensivity of our predi ction s with respect to the pa- 
rameterization (119) and (|l22j ) we also plot in Fig. || our 
results for the same cosmologies when we use K — 0.88 
(this would correspond to an initial power spectrum in- 
dex n = —1.5) in (122). We can see in the figure that 
our predictions are not too sensitive to k (for reasonable 
values of k) within the range —0.015 < M ap < 0.015 (the 
difference would look even smaller in Fig. ||). In partic- 
ular, the variation with k of our results is much smaller 
than the difference between both cosmologies. This sug- 
gests that one could use the deviation of the PDF with 
respect to a Gaussian to estimate the cosmological pa- 
rameters. A well-known tool to measure this signature is 
the skewness but one could devise other statistics which 
would take advantage of the expected shape of the PDF to 
maximize their dependence on cosmology. However, such 
a study is beyond the scope of this article. We can see 
in Tab.|2| that we underestimate somewhat the skewness 
of P(M ap ). However, it is not clear whether this is due 



to the parameterization ( |119| ) or to the use of HEPT in 
(122). We note that the tail of P(M ap ) for large negative 



values of M ap appears to be slightly more sensitive to n 
than the tail at positive M ap . This could be related to the 
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Table 2. Cosmological models and results obtained with Us filter. The sixth and seventh lines show the variance 
and the skewness of P(M ap ) from numerical simulations (Jain, Seljak & White 2000). Results for sCDM and ACDM 
models have been directly taken from Reblinsky et al. (1999) and have been obtained with one realization only. The 
eight and ninth lines show the reconstructed PDF properties. 





sCDM 


rCDM (5 realizations) 


OCDM (7 realizations) 


ACDM 






1 


1 


0.3 


0.3 




Oa 











0.7 


H [km/s/Mpc] 


50 


50 


70 


70 






0.6 


0.6 


0.85 


0.9 




r 


0.5 


0.21 


0.21 


0.21 


° M a p , 


9 = 4' 


0.00730 


0.00542 ± 0.00017 


0.00449 ± 0.00017 


0.00495 


s?»°, 


9 = 4' 


159 


258. ± 38. 


545. ± 30. 


347 


° M a p , 


9 = 4' 


0.00721 


0.00497 


0.00469 


0.00477 




9 = 4' 


148 


180 


395 


300 



fact that the behaviour of P(M ap ) for Af ap < — crjvf ap i s 
more sensitive to the detailed properties of the p— point 
correlation functions (see Valageas 2000c for a study of 
this point). 

Finally, we note that although 9 — 4' corresponds to 
non-linear scales it is not very far from the quasi-linear 
regime. However, the aperture mass M ap probes the non- 
linear density field for filters with larger angular scale 
than the convergence k. Indeed, since the aperture mass 
involves compensated filters the contribution from low-fc 
modes is strongly suppressed (see Fig. so that P(M ap ) 
is governed by the properties of the density field at the co- 
moving wavenumber k ~ 2/(T>9). In contrast, the conver- 
gence k shows a more important contribution from larger 
wavelengths which implies that in order to probe non- 
linear scales only one must set the filter size 9 farther away 
into the small-scale non-linear regime. This also means 
that in principle the aperture mass could be a more con- 
venient tool than the convergence since it should be easier 
to separate the non- linear and quasi- linear regimes, while 
for an important range of angular scales the convergence 
should be sensitive to the transitory regime between both 
domains. However, a possible caveat is that the statistics 
of the aperture mass depend on the detailed behaviour of 
the p— point correlation functions (and not on their av- 
erage over spherical cells only), and therefore requires a 
better understanding of them. 

As discussed above, this property of the aperture mass 
to probe a narrow range of wavenumbers makes it easier 
to avoid the intermediate regime as one can select obser- 
vation windows which arc cither in the quasi-linear or in 
the strongly non-linear regime. However, it would clearly 
be interesting to obtain a model which would also cover 
this transitory range. Unfortunately, this is rather difficult 
as one cannot use the simplifications which appear in the 
two extreme regimes. An alternative to a rigorous calcu- 
lation would be to use an ad-hoc parameterization which 
would smoothly join the quasi-linear regime to the highly 
non-linear regime, for instance in the spirit of HEPT as 



described in van Waerbeke et al. (2000b) for the skewness 
of the convergence (note that our model for non-linear 
scales is based on a simple ansatz which is not rigorously 
derived). However, although the quasi- linear and strongly 
non-linear regimes share the same gross features, like the 
relation ([34j) which describes the projection effects, their 
detailed properties are different. In particular, although in 
both cases we have the scalings (|3^) the correlation func- 
tions obtained in the quasi-linear regime are not given by a 
tree-model as in ( |100| ). This leads to the difference between 
the two- variable system (|77|)-(|80|) and the integral relation 
(117)- (118). Neverthele ss, a s imp le prescription would be 
to recast the relations (|il^)-(|iTgl) into the form (|77|)-(|o|) 
by approximating the integrals over ■& by the difference 
between two mean values, corresponding to the inner and 
outer regions of the filter Ubv and characterized by two 
averages t± and ti - Then, the shift from the quasi-linear to 
the highly non-linear regime would simply be described by 
a smooth interpolation of the generating function C( T )j i-e. 
of the sole parameter n. However, such a study is left for 
a future work as in this article we prefered to lay out the 
formalism needed to study the statistics of the aperture 
mass and to focus on the two regimes which have already 
been tested in details against numerical simulations, so as 
not to introduce a new specific parameterization. 



7. Conclusion 

In this article we have described methods that allow exact 
reconstructions of the one-point PDF of the local aperture 
mass in weak lensing maps. These methods do take into 
account the projection effects but not all the nonlinear 
couplings between the local density field and the observed 
distortions field such as lens-lens coupling effects, or de- 
parture from the Born approximation. 

In the course of this paper we have examined both 
the quasi-linear and non-linear regimes. In particular, al- 
though the details of the calculations are specific to each 
case we have pointed out the generic properties which are 
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Fig. 8. The aperture mass PDF for a source at redshift z s = 1 for 4 cosmologies and with the angular window 9 = 4'. 
The dashed line shows the Gaussian which has the same variance. The points show the results of N-body simulations 
from Jain, Seljak & White (2000). Results for sCDM and ACDM models have been directly taken from Reblinsky et 
al. (1999). 



common to both regimes and the features brought about 
by the projection effects. For instance, in both quasi-linear 
and non-linear domains the PDF P(M ap ) should show two 
asymmetric exponential tails. Our methods are quite gen- 
eral and can be extended in a straightforward fashion to 
other statistics. In the quasi-linear regime our approach 
can be applied to any filter which is axisymmetric while 
in the non-linear regime there are no restrictions. In par- 
ticular, in this latter case our results can be extended to 
multivariate statistics (which can be obtained from filters 
which consist of several disconnected parts). 

We have briefly investigated the dependence of the 
PDF P(M ap ) with the shape of the filter. Thus, we have 
checked that for filters with a large compensation radius 
we recover approximatly the shape of the PDF P(k) which 
is relevant for the top-hat filtered convergence. 

Finally, we have checked that our predictions agree rea- 
sonably well with the results of available numerical sim- 
ulations (although we have not included any noise effect 
at this level) at scale about 4' where the data should pro- 



vide the largest signal to noise ratio (e.g. Jain & Seljak 
1997). In particular, we recover the asymmetric shape of 
the PDF. Moreover, our approach provides a prediction 
for the full shape of the P(M ap ) while earlier models were 
restricted to the positive tail of the PDF. We have also 
shown that the difference between the PDFs obtained for 
two cosmologies (tCDM and OCDM) is larger than the 
inaccuracy of our predictions (due to parameterization we 
need to introduce to describe the underlying 3D density 
field). This suggests that our results could be used to es- 
timate the cosmological parameters. Thus, in addition to 
the skewness which is traditionally used to this purpose 
one could take advantage of the expected shape of the 
PDF to build other statistics which would maximize the 
dependence on the seeked parameters. Such a study is left 
for further work. 
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Appendix A: Order of the singularity of tpproj.Cy) 

As noticed in Valageas (2000a) the exponent ui StC of the sin- 
gularity of ip cy i,(y) (as defined in (^)) translates into the ex- 
ponent u! StP — ui StC — 1/2 for the projected generating function 
(/Cproj. (j/)- However, in the cases encountered in this article we 
have uj s .c = —3/2 for both the quasi-linear and highly non- 
linear regimes. Then u S:P = —2 is an integer but the generating 
function ip P roj.(y) is still singular at the points y£!jf J through 
logarithmic factors. To see this, it is convenient to take the 
third derivative of the relation (|34|) which is governed by the 
singularity at y^™ 1 and diverges for y —> y^™ 1 (while the lower 
derivatives of V'proj^J/) remain finite at y^° J ). This yields: 



Coj.(y) = j dTiF(iiyMKyrt;>.[yF(n)MK)]. (A.i) 

For y — > y^™ 1 the integral is dominated by the values of 1Z 
around the point where the factor F(lZ)i})e{lZ) is maximum, 
since tpty di verg es as \y — y^^l^ 3 ^ 2 at this point. Thus, we 
obtain from (A.l): 



V^vlt ■■ ^j.(v)~ / ATI \(l-n 2 )y-y™\ 3/2 .(A.2) 



After the change of variable t = \y/(y — yg± y )\7l 2 we obta 



proj (3) / \ f°° dt . jA-3/2 M 1/2 

Jo V* IV-Vsil 

which gives: 



\y-yT± } \ 



Finally, the integration of this relation leads to: 

y -* ylT ■ VproiXv) ~(y- y p s ±') 2 - vT± J l 

where we only wrote the most singular term. 



(A.3) 



(A.4) 



(A.5) 



